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We analyze the photovoltaic current through a double quantum dot system coupled to a high- 
quahty driven microwave resonator. The conversion of photons in the resonator to electronic excita- 
tions produces a current flow even at zero bias across the leads of the double quantum dot system. 
We demonstrate that due to the quantum nature of the electromagnetic field in the resonator, the 
photovoltaic current exhibits a double peak dependence on the frequency uj of an external microwave 
source. The distance between the peaks is determined by the strength of interaction between photons 
in the resonator and electrons in the double quantum dot. The double peak structure disappears 
as strengths of relaxation processes increases, recovering a simple classical condition for maximal 
current when the microwave frequency is equal to the resonator frequency. 

PACS numbers: 73.23.-b, 73.63.Kv, 42.50.Hz, 73.50.Pz 



I. INTRODUCTION 

The interaction of electrons in conductors with electro- 
magnetic fields has long been considered within a classi- 
cal picture of electromagnetic (EM) radiation. A widely- 
known example is the photon assisted tunneling (PAT) 
in double quantum dot (DQD) systems^ when the EM 
field brings an electron trapped at the ground state to an 
excited state and facilitates electron transfer. This clas- 
sical description of the EM field breaks in high-quality 
microwave resonators based on superconducting trans- 
mission line geometry.— Interaction of such EM fields with 
electronic devices require a quantum treatment known as 
the circuit quantum electrodynamics (cQED)»^— 

Recently, several experimental groups studied systems 
consisting of a superconducting high quality resonator 
and a DQD^"— or a voltage biased Cooper pair box.— The 
coupling strength between a resonator photon mode and 
electron states in a DQD is characterized by the vacuum 
Rabi frequency g with reported values in the range of 
g/2Tr ^ lO^Hz. These systems call for re-examination of 
the PAT by taking into account a quantum description of 
the EM field in terms of photon excitations. One may ex- 
pect at least two important distinctions from the classical 
treatment: (1) the Lamb shift that renormalizes quantum 
states of electrons and photons; (2) spontaneous photon 
emission that breaks symmetry between absorption and 
emission processes and is important in systems with ei- 
ther a finite voltage bias between the leadsii"— or an 
inhomogeneous temperature distribution^^ 

In this paper we study the photovoltaic current 
through a DQD coupled to a high-quality microwave res- 
onator at zero bias across the DQD. The resonator is 
driven by external microwave source that populates a 
photon mode of the resonator, sec Fig. [TJ The pho- 
tons excite electrons in the DQD and produce electric 
current even at zero bias, similar to the classical PAT 
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We show that due to the coupling of electrons 
and photons, the current as a function of the source fre- 
quency has a multiple peak structure with splitting be- 
tween the peaks determined by the coupling strength g 



//w resonator 



1^ 



(b) Z 





























e 






no. 


1 




























\9) 







BL PL BM PR BR 



FIG. 1: (Color online) (a) An illustration of a DQD and a 
transmission line resonator coupled to an external microwave 
source ^w. (b) A schematic view of the DQD. Electrons are 
confined to the left (L) and right (R) dots by barrier gates 
BL, BM, and BR that also control electron tunneling rates 
between the source, S, and the left dot, the left and right dots, 
and the right dot and the drain, D, respectively. Electrostatic 
energies of two quantum dots are defined by the plunger gates, 
PL and PR, and the PL gate is also connected to an antinode 
of the resonator, see e.g. Refs.QB (c) Electronic states of the 
DQD are presented in both the eigenstate basis (solid lines) 
and the left -right basis (dashed lines). Tunneling from the 
excited state, |e), to the left/right lead, with rate T^^i^/ji and 
from the left /right lead to the ground state, \g), with rate 



L/R,g 



are illustrated by arrows. 



and reflects the Lamb shift of electronic energy states. 
We also demonstrate that the interaction-induced split- 
ting is sensitive to the energy and phase relaxation rates 
in the DQD. 

We note that the photovoltaic effect discussed here is a 
common phenomenon when the current in an electronic 
circuit is generated by out-of-equilibrium EM environ- 
ment. Examples of this phenomenon include the current 
response of a DQD in the vicinity of the biased quan- 
tum point contaclj^ or another circuit element out of 
equilibrium^^ with the electronic system. However, be- 
cause out-of-equilibrium photons of the environment have 
a broad spectrum, the generated current does not exhibit 
a resonant dependence on parameters of the system that 
we observe in a system of a single mode high quality res- 
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onator and a DQD. 

II. MODEL 

We consider a DQD system with each dot connected 
to its individual electron reservoir at zero temperature 
and at zero bias between the reservoirs, see Fig. [Hb,c). 
The gate voltages of the DQD are adjusted near a triple 
point of its stability diagramii To be specific, we choose 
a triple degeneracy point between (TV;, Nr), {Ni + 1, Nr) 
and {Ni, Nr + I) electron states in the DQD and denote 
these states as |0), \L) and \R), respectively. We model 
the system by the Hamiltonian H = iJoQD + Hr + ^^int , 
where Hdqd describes states with an extra electron in 
the left or right dot, \L) or \R): 

Hbqd = ^£Tz + Tt^, (1) 

with e being the electrostatic energy difference between 
the two states, and T being the tunnel matrix element 
of an electron between the dots. The Pauli matrices are 
defined in the subspace of states \L) and \R) as ~ 
\R) {L\ + \L) {R\ and = \R) {R\ - \L) {L\. A resonator 
driven by an external microwave source is described by 
the Hamiltonian 

Hr ~ hbJoa^ a + 2hF{a^ + a) cos uit (2) 

with a (a^) denoting the annihilation (creation) opera- 
tors for microwave photons in the resonator, HF being 
the amplitude of the external drive of the resonator and 
ujQ (cj) being frequency of the resonator (source). The 
interaction between the microwave field and the DQD 
system is represented hy^ 

Hint = fi9o{a^ + a)T^. (3) 

This interaction describes the shift of energy difference 
between states \R) and \L) due to the electric potential of 
the plunger gates defined by the microwave photon field. 
We assume that the photon field is distributed between 
the left and right plunger gates, see Fig. [Ijb) and does 
not influence the source and drain voltage to avoid the 
rectification effects.—"— 

Further calculations are more convenient in the basis 
of the ground, \g), and excited, |e) states of the Hamil- 
tonian, Eq.ll]): 

|e) = cos(6l/2) \L) + sin(6'/2) \R) , 

|5) = -sin(0/2)|L)+cos(0/2)|i?). (4) 

Here = arctan(2T/e) characterizes the hybridization of 
the \L) or \R) states. The energy splitting between the 
eigenstates Ml = Ve^ + ^T^ can be tuned independently 
by varying e and 7" via dc gate voltages. We further elim- 
inate the time-dependence in Hamiltonian Eq.Q by ap- 
plying unitary operator U = exp{—iujt{a''a + (Tz/2)) and 



utilize the rotating frame approximation to obtai n^ i^ 

+ (cjQ — (jj)a^a + (7(a(T^ + a^a ) + F{a^ + a), 

where g = go sin 9 characterizes the actual strength of the 
coupling between the microwave field and DQD states 
responsible for photon absorption or emission, the Pauli 
matrices CTz ~ \e) (e| — \g) {g\, = \g) (e| and a'^ = 
|e) {g\ are defined in terms of eigenstates of the electron 
Hamiltonian, Eq. ([T]). 

We analyze the behavior of the system with Hamil- 
tonian Eq. ([5]) in the presence of relaxation in electron 
and photon degrees of freedom by employing the Born- 
Markov master equation for the full density matrix 

p = CtotP = --[H,p]+ CdissP- (6) 

The first term on the r.h.s. of Eq.® describes the uni- 
tary evolution of the system and the second term ac- 
counts for the dissipative processes in the resonator and 
DQD systems^ 

CdissP = n'D{a)p + jT>{a~)p+'^'D{<T^)p 

2 (7) 

+ {TL.g + TR^g)V{cl)p + (Fe.L + Fe,fl)I?(Ce)p, 

where 'D{x)p = (2xpx^ — x^xp — px'^x) /2 is the Lind- 
blad superoperator. The relaxation of the photon field 
in the resonator with rate k is represented by K'D{a)p and 
the electron relaxation from the excited state |e) to the 
ground state \g) with rate 7 is represented by ^'D{a~)p. 
The last two Lindblad superoperators account for the 
loading of the ground state \g) and unloading of the ex- 
cited state |e) of the double quantum dot via electron tun- 
neling in terms of operators Ce = |0) (e| and = \g) (0|, 
respectively. The tunneling rates F^^g = F; cos^(6'/2), 
Ffl,g = F^sin2(6l/2), T^x = Tisin^{e/2) and Te,R = 
Tr cos^(0/2) arc written in terms of tunneling rates F;/^ 
in the basis of \L) and \R) states. 

Note that in Eq.®, the dynamics of state |0) only 
appears via the tunneling terms involving D{ce)p and 
T>{c^g)p. These terms can be categorized by whether the 
empty state is loaded from the left or right lead with co- 
efficients depending on projection of the eigenstates onto 
the left/right states, as shown in Fig. [T] In this picturep^ 
the photovoltaic current is given by 

I = eTr (^cos2 ^ (e| p |e) - sin^ ^ (0| p |0)^ . (8) 

in terms of the reduced density matrix p = Trph{p}, 
where we traced out photon degrees of freedom of the 
resonator. We also analyze the number of photons in the 
resonator, 

N = Tr{a^ap}, (9) 

where we trace out both photon and electron degrees of 
freedom. 
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FIG. 2: (Color online) The average number of photons in 
the resonator and the photovoltaic current as functions of 
level bias e for T/2-k = 1 GHz, F = 50 ^s~^ and ujo/2-k — 
8 GHz. The current is generated near the resonant condition 
when e — zt^JfALu"^ — 4T^ (vertical lines). The three curves 
represent different dephasing rates 7^ the DQD. 



III. RESULTS 

The average number of photons in the resonator, iV, 
and the dc component of photocurrent can be found using 
the steady state solution of the master equation, ^ , with 
p = 0. We numerically find the full density matrix p for a 
double quantum dot and photon field of the resonator in 
the Fock's space using Quantum Optics Toolbo5i22, and 
QuTiP— , both of which provide consistent results. The 
steady state solution for the density matrix p defines the 
average number of photons N, Eq. ([5]), and the pho- 
tocurrent, Eq.®. 

Our choice of parameters is motivated by Ref . . We 
choose the relaxation rate 7 = 27r x 25 MHz, the res- 
onator relaxation rate K/27r = 8 MHz, tunneling ampli- 
tude between the individual dots T /2tt = 1 GHz, the 
tunneling rate from a dot to a lead Vi/r = 27r x 30 MHz, 
and the resonator frequency ujqI^-k = 8GHz. We note 
that to keep the coupling constant finite, we have to take 
T ~ ft^l, since g = go sin 6*, Eq.®, vanishes for T = 0. 
Below we fix go/27r 200MHz. 

First, we investigate dependence of the photocurrent 



FIG. 3: (Color online) The average number of photons in 
the resonator and the photovoltaic current as a function of 
the frequency tj of the microwave drive for T/2ti — 1 GHz, 
F = 50 fis'^ and n = cjo = 27r X 8 GHz, g/2Tv = 48.5 MHz. 
For 7^ = 0, both average number of photons A*' and the photo- 
voltaic current show local minima at cj = a;o and local maxima 
near oj = (^-i.i — Eo)/h, shown by vertical lines. As the de- 
phasing rate 7^ increases, the double peaks merge to a single 
peak at lj = ujo. 



on the separation between energy levels in the double 
quantum dot, controlled by the electrostatic energy dif- 
ference e. We take frequency uj of microwave source to 
be equal to the resonator frequency, a; = wq, and fix the 
drive amplitude F = 50 p,s~^. Dependence of the average 
number of photons in the resonator and the photocurrent 
on energy e is presented in Fig. [2] for three values of the 
dephasing rate j^/2tt = 0, 10, 20MHz. As the energy 
difference between the excited and ground states of the 
quantum dot goes through the resonance Q ~ ujq, we ob- 
serve a significant suppression of the photon number in 
the resonator, see the top panel and the inset in Fig. [S] 
This is expected behavior because the DQD system en- 
hances photon absorption in the resonator at f2 — wo- 
Absorbed photons cause transitions between the ground 
and excited electronic states. These electrons tunnel to 
the leads and generate electric current though the DQD. 
This current is shown in the lower panel of Fig. [5] and is 
peaked at e = ±^/fK^^~4n^ or e/{2TTh) ~ ±7.75GHz, 
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FIG. 4: (Color online) The average number of photons in the 
resonator and the photovoltaic current as a function of the 
frequency uj of the microwave drive for detuned DQD and 
resonator system with n/27r = 8.1 GHz, a;o/27r = 8.0 GHz, 
the intradot tunneling T/2n = 1 GHz, and the drive ampli- 
tude -F = 50 MHz. The photon average number has a peak at 
UJ = (-El,- — Eo)/h, Eq. (|10p . while the photovoltaic current 
exhibits a double peak feature at w = (i5i,± — Eo)/h (vertical 
lines). 



indicated by dashed vertical lines. 

One feature in Fig. [2] is that the photon number is also 
reduced at zero bias e, when the photovoltaic current 
vanishes. This suppression is a result of strong enhance- 
ment of the coupling constant g = go at e = 0, resulting 
in stronger dissipation in the resonator and increase of 
off-resonant absorption rate. At the same time the pho- 
tovoltaic current vanishes at e = due to cancellation 
between the two terms in Eq.®. 

The curves for the photon number and the current do 
not significantly change after the dcphasing rate 7^ is 
introduced in addition to the energy relaxation rate 7. 
Dephasing smears the resonant condition for the pho- 
ton absorption by the DQD and has two effects: (1) the 
number of photons increases a little near the resonance 
ri ~ Wo, see the inset in Fig. [21 (2) the resonant absorp- 
tion of photons by electrons is suppressed resulting in 
reduction of the photocurrent. We note that in the case 
presented in Fig. [2] the first effect is stronger than the 



second effect and dephasing increases the magnitude of 
photocurrent for the case of fixed w = wq. 

Next, we consider the case when the frequency of the 
microwave source, w, is varied while the energy splitting 
HQ of the DQD and the resonator frequency cjq are fixed. 
The microwave radiation is mostly reflected when its fre- 
quency does not match the difference between energies 
£"„ _!_ of the resonator and DQD system defined by the 
Jaynes-Cummings spectrum: 



En 



(10) 



where A = wq — is the detuning between the DQD and 
the resonator. We demonstrate that for DQD with weak 
energy and phase relaxations, this resonant admittance 
of the microwave source to the resonator results in the 
peak structure of the average photon number and the 
photocurrent. 

In Fig. [21 we plot the average number of photons in 
the resonator and the photocurrent as a function of the 
drive frequency to for ujq = ^ and for the choice of other 
system parameters identical to those for curves in Fig. [51 
At vanishing dephasing rate, 7^ = 0, we observe a double 
peak feature in both photon number and photocurrent 
curves, see Fig. [3 These peaks at uj± = {Ei,± — Eo)/h 
are defined by the level spacing of the Jaynes-Cummings 
Hamiltonian and are shown by vertical dashed lines in 
Fig. [31 The two peaks merge at w = wq as the dephas- 
ing rate increases and destroys quantum entanglement 
between photons and DQD states. 

At finite detuning between the resonator and the DQD, 



A > 



27r X 48.5 MHz, the eigenstates of the system be- 



come dominantly photon states or electron states of the 
DQD. As a result, the microwave source increases the 
number of photon excitations in the resonator when the 
microwave frequency is in resonance with the transition 
between the photon-like states, = (i?i,- — En)/h. 

But the source has a weak effect at the resonance with the 
electron-like states at frequency wi^-i- = {Ei^+ — Eq)/H. 
We present the corresponding dependence of the pho- 
ton number and the photocurrent in Fig. 31 for a;o/27r ~ 
8 GHz, n/2TT = 8.1 GHz (A = 100 MHz) and other pa- 
rameters identical to those for in Figs. [Hand [31 We in- 
deed observe one large peak in the photon number near 
the resonant condition for the dominantly photon state 
with energy while the photon number does not show 
significant enhancement near the second resonance, cor- 
responding to the transition to the dominantly electronic 
state with energy £'1,+. The photocurrent still exhibits 
double peak feature, but the peak corresponding to the 
photon resonance is higher, when the microwave drive 
produces a higher photon population. 

We now consider a more idealistic regime of signifi- 
cantly reduced tunneling and relaxation rates F; = F,. = 
7 = 27r X 100 kHz, the drive amplitude F/27r = 30 MHz 
and Wo = = 27r X 8 GHz. In this case additional reso- 
nances develop, see Fig. [SI These resonances correspond 
to excitations of several photons in the cavity by the mi- 



5 



I 0.2 

3 

o 



^ 15 



a 

3 
U 

73 J 



(a) A 

- — — 1 — 


1 ' 1 

> 


L 


J 


1 1 


I 



-0.1 

(b) ' 

0.5 




-0.05 





)/2: 
(C) 



0.05 



0.1 



{(i)-&^)/2Ti, GHz 



Jl 



12 3 

n 




IO,g> 



FIG. 5: (Color online) (a) The average number of photons 
in the resonator and the photovoltaic current as a function 
of the frequency iu of the microwave drive for Q — ujo ~ 
2n X 8 GHz, the intradot tunneling T/2tt — 1 GHz, and the 
drive amplitude F/2tt — 30 MHz and extremely low tunneling 
rates to the leads and the energy relaxation rate, F; = Fr = 
7 = 27r X 100 kHz. The photon average number and the 
photocurrent have several peaks at ujn,± = {En,± — Eo)/nh 
with n = 1, 2, 3, these frequencies, calculated from Eq. pop 
are shown by vertical lines), (b) The histogram presents the 
probabilities Pn to have n photons in the resonator steady 
state at drive frequency uji (dark bars) and 102 (grey bars) . (c) 
The diagram is a schematic picture for the Jaynes-Cummings 
energy levels showing single and two photon excitations. 



crowavc source. When the frequency of the source satis- 
fies hujn = En^± — Eq, the DQD-resonator system expe- 
riences transitions from the ground state to the energy 
state i?„.±, cf. Ref. [l^- These multiphoton transitions 
result in peaks of the average photon number and the 
magnitude of the photocurrent. Curves in Fig. O have 
three pairs of peaks at frequencies ujn.± = '^o i g/\/n 
marked by vertical dashed lines for n = 1,2,3. We no- 
tice that for CO = U!i^2 the average photon number is 
nearly the same, see the top panel in Fig. El^a), while 
the photon distribution function is different. Fig. Hfb): 
at uj = UJ2,- a non-zero P2 develops for a probability 
that the resonator contains two photons. This difference 
in P„ indicates that the microwave drive line does not 
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FIG. 6: (Color online) Dependence of the average photon 
number and the current on microwave frequency a; for several 
values of the drive amplitude F/2n = 5, 10, 15 MHz, at zero 
dephasing 7,^ = and other parameters are the same as in 
data in Fig. [S] As the amplitude of the drive increases, the 
two peaks merge together to a single peak at lu = ojo- 



match the resonator to produce a two photon occupation 
of the resonator at a; = wi , but it matches the resonator 
to populate the state with the energy i?2,±, which then 
decays to the lower energy states with 71= 1,0. 

Next, we investigate dependence of the photon number 
in the resonator and the magnitude of the photovoltaic 
current for different amplitudes F of the drive. The above 
discussion was mostly focused on a resonator containing 
less than one photon. As the drive increases, the double 
peak feature evolves to a single peak at the drive fre- 
quency equal to the frequency of the resonator, lo = ujq. 
We interpret this cross-over as a signature of changed 
hierarchy of the terms in the system Hamiltonian. At 
weak drive, we have a JC Hamiltonian with its peculiar 
energy levels, Eq. pH)) . and the drive can be viewed as a 
weak probe testing the spectral structure of the coupled 
resonator and DQD system. Once the drive reaches the 
strength of the g coupling, g ~ 27r x 50 MHz, a proper 
way to treat the system is to start with the Floquet-type 
state a^^'^^1^^ of the driven resonator and then to take into 
account the interaction of these states with the DQD sys- 
tem as a perturbation. In this picture, the photon reso- 
nance happens at a; = wq. The coupling g is responsible 



6 



for the formation of the broader "wings" in curves for 
the average photon number and the photocurrent. These 
wings are more pronounced in the photovoltaic current, 
which is entirely due to the coupling between resonator 
and DQD. This broad structure of the generated current 
as a function of the source frequency is preserved even at 
stronger drive. Thus, the shape of the photovoltaic curve 
might provide an experimental approach to quantify the 
strength of the JC coupling constant. 

IV. DISCUSSION AND CONCLUSION 

We analyzed the photovoltaic current through a DQD 
system at zero voltage bias between the leads. The dou- 
ble quantum dot interacts through its dipole moment to 
a quantized electromagnetic field of a high quality mi- 
crowave resonator. The interaction is described by the 
Jaynes-Cummings Hamiltonian of a quantized electro- 
magnetic field and a two level quantum system, repre- 
sented by ground and excited electronic states of the 
double quantum dot. When a weak microwave radiation 
is applied to the resonator, the source acts as a spec- 
tral probe that causes excitation of the system when the 
energy difference between its eigenstates is equal to the 
photon energy hw of the source. If this resonance con- 
dition is satisfied, the microwave source populates the 
photon mode of the resonator and generates a direct cur- 
rent though the double dot system even at zero bias. 

We demonstrated that at finite, but still low energy 
and phase relaxation rates of the DQD, both the average 
number of photons in the resonator and the photocur- 
rent through the DQD have a double-peak structure as 
functions of the frequency of the microwave source. This 
double peak structure reflects an avoided crossing of the 
energy states of the DQD and the resonator photons due 
to the interaction between the two subsystems and is 



reminiscent of the Lamb shift by a single electromagnetic 
mode. We also found that in the limit if extremely weak 
relaxation rates of the DQD, multiphoton resonances de- 
velop when the energy difference between the states of 
the coupled system is a multiple of hw. 

As energy and phase relaxation rates of the DQD in- 
crease, the peaks in the photon number and the photocur- 
rent broaden and eventually merge in a single resonance 
peak at the frequency wo of the resonator. In this limit, 
the resonator photon mode and the DQD arc no longer 
described as an entangled quantum system and the reso- 
nant condition for the interaction of the microwave source 
with the system corresponds to equal frequencies of the 
source and the resonator mode, u ~ ujq. 

At stronger microwave drive, frequency dependence of 
the average photon number in the resonator evolve from 
the Jaynes-Cummings double peaks aX uj = ujo ± g to & 
single peak at the resonator frequency wq. The single 
peak at w = Wo is a result of multi-photon transitions 
at strong drive by the microwave source that all merge 
together due to finite width of multi-photon resonances. 
Similar evolution to a single peak occurs for the pho- 
tocurrent response, although the photocurrent curve has 
a broader width as a function of the source frequency w, 
this width corresponds to the strength of the coupling g 
between the photon mode of the resonator and the DQD 
and may be used to characterize the strength of this cou- 
pling in experiments. 
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